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Abstract 

Half-filled two-leg Hubbard ladders have spin-gapped short-range anti- 
ferromagnetic correlations while three-leg ladders have power law antofer- 
romagnetic correlations, and both systems have d 2 .2_ J/ 2-power law pairing 
correlations when they are doped. Thus these ladders exhibit some of the 
phenomenology seen in the layered cuprates. Here we report results for half- 
filled frustrated Hubbard ladders, based upon ladder segments taken from a 
tight-binding model of k-BEDT-TTF. Although these ladders are half-filled, 
varying the degree of frustration can drive them across an insulator-metal 
transition. We suggest that the spin, charge and pairing correlations of these 
frustrated ladders near the insulator-metal transition provide support for the 
notion that k-BEDT-TTF is a strongly correlated superconductor. 

The organic k-BEDT-TTF crystal consists of weakly coupled two-dimensional layers 
of BEDT-TTF molecules separated by insulating layers of anions such as Cu[N(CN)2]Br. 
Under pressure, and at low temperature (T ~ 10K), this system undergoes a phase tran- 
sition from an insulating antiferromagnetic phase to a superconducting phase. Based upon 



parameters obtained from quantum chemistry calculations, it is believed that this material 
is a strongly correlated electron system. Within this framework, one can seek to under- 
stand the physics of the low temperature k-BEDT-TTF phase diagram in terms of a model 
containing a strong on-site Coulomb interaction and a tight binding parameterization of the 
relevant band structure. 

A diagram showing the arrangement of the k-BEDT-TTF molecules and their dominant 
hopping integrals, as proposed by Kino and Fukuyama |§, is shown in Fig. 1. In this model 
each lattice site, denoted by a circle, contains two BEDT-TTF molecules which are strongly 
hybridized leading to a large splitting between their bonding and antibonding orbitals. The 
three electrons on each site leave the antibonding orbital with one electron and it is the 
resulting half-filled t — r antibonding band with a strong on-site Coulomb interaction U that 
has been proposed as a minimal model. For r = 0, this is just a near neighbor 2D Hubbard 
model which at half-filling has an insulating antiferromagnetic ground state. When U is 
large compared to the bandwidth, and r increases, the diagonal exchange leads to frustration 
of the antiferromagnetic order and the possibility of a spin gapped phase. Alternatively, for 
smaller values of U, the system is expected to go from an insulating antiferromagnetic phase 
to a paramagnetic metal and possibly superconducting phase as r increases. KJ 

Motivated by this scenario and the geometry of the tight binding lattice, we report results 
obtained from density matrix renormalization group || (DMRG) calculations on the two 
different ladders indicated by the large arrows shown in Fig. 1. The horizontal ladder, 
corresponding to a one-dimensional chain with a near neighbor hopping t and a next near 
neighbor hopping r has been previously studied || and we will summarize its properties in 
our conclusion. 

The diagonal ladder of Fig. 1 has a Hamiltonian 

H = -t Y ( c l*Cj,« + h -c) - r Y ( c l<r c j,a + h.c) + UY n i\ n H (!) 

<ij>,o [ij], a i 

where cj CT (cj )(J ) creates (annihilates) an electron on site i with spin a. The sum < ij > 
runs over pairs of nearest neighbor sites, and the sum [ij] runs over pairs of lattice along 
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only one diagonal of the ladder. In the following we measure energies in units of t, and the 
diagonal hopping r is positive. To investigate this system we applied the density-matrix 
renormalization group Q technique for system up to size L y x L x = 2 x 32 with open 
boundary conditions. We kept up to 800 states so that the maximum discarded weight was 
lO" 6 . 

The charge gap 

A c = l - [E (N + 2,S = 0) + E (N - 2, S = 0) - 2E (N, S = 0)] (2) 
and the spin gap 

A s = E Q (N, S = l)- E (N, S = 0) (3) 

for a half-filled (N = 64) 2 x 32 ladder are shown in Figs. 2a and 2b, respectively. Here 
Eq(N, S) is the ground-state energy for a state with N electrons and spin S. The gaps are 
plotted versus U for different values of r. They were determined by finite-size scaling of the 
result from systems of size L x = 8, 16, 24 and 32. The charge gap is set by some fraction 
of U and hence is relatively easy to determine. The spin gap, which is set by the effective 
exchange interaction, is small even for the unfrustrated r = ladder. When r is finite, the 
resulting frustration further reduces the spin gap. Hence for small values of U, there are 
strong finite size effects which make it difficult to determine the spin gap very precisely. But 
here we are interested in the qualitative phase diagram rather than quantitative. We have 
also calculated the charge and spin gaps for r = 1.2 and 1.5, which were both finite for small 
U. 

For r = 0, we have a half-filled 2-leg Hubbard ladder and as expected both the charge 
and spin gaps open up as U turns on. In opposite, for 0.7 ~ r ~ 1.1, a finite value of U 
is required to open up a charge gap. For example, for r = 1, a value of U ~ 4 is required 
before the system develops a charge gap and becomes insulating. When r is present, it leads 
to a frustrating diagonal exchange interaction. For a sufficiently large value of r (r ~ 1), 
the dominant singlets are formed along the r diagonal bond rather than the t leg and rung 
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bonds of the ladder. Based upon the results of Fig. 2, a schematic phase diagram for the 
"diagonal" ladder is shown in Fig. 3a. 

The Hamiltonian of the "horizontal" ladder shown in Fig. 1 is identical to that of the 
"diagonal" ladder for r — 1. However, for r^l, the Hamiltonians for the two ladder differ. 
In particular, for r = 0, the "horizontal" ladder reduces to the 1-leg Hubbard model which at 
half-filling has a charge gap but no spin gap as opposed to the diagonal ladder which reduces 
to the 2-leg Hubbard ladder with a spin gap. Thus for small values of r, the "horizontal" 
ladder is characterized by power law antiferromagntic correlations and a charge gap. The 
phase diagram for the "horizontal" ladder obtain in Ref. H, is shown in Fig. 3b. In this case 
one has an insulating phase with power law antiferromagnetic correlations at smaller values 
of r and a spin gapped insulator with short range antiferromagnetic correlations separated 
from a smaller U metallic phase at larger values of r. 

In order to explore the tendency for pairing, we have measured the rung-rung singlet 
pair field correlation function 

D{i) = (A l+e Aj). (4) 

The operator 

A» = 4l,T C i2,| — C il,l C i2,1 (5) 

creates a singlet pair on the i rung and A i+ £ destroys it on the (i + £) rung. In Fig. 4 we 
show results for D(£) versus £ for r = 1 and values of U above and below the insulator-metal 
transition. For U = 6, the system is in the insulating phase and D(£) decays exponentially. 
When U decreases below 4, the charge gap vanishes and the pair field correlations exhibit a 
power law decay. 

In summary, DMRG calculations for 2-leg ladder segments of the model k-BEDT-TTF 
lattice give phase diagrams which have features similar to the proposed phase diagram of 
the 2D system. |TJ In particular, for small values of r, they both have an insulating phase 
with antiferromagnetic correlations. In the diagonal ladder, these correlations are short 



range and the system is spin gapped while for the horizontal ladder these antiferromagnetic 
correlations can be power law or short range depending upon the size of r. In addition for 
moderate to smaller values of U, both ladders can go into a metallic state as r increases and 
this state exhibits power law singlet pairing correlations. Thus one could imagine that as 
the pressure is increased the system crosses from an insulating state with antiferromagnetic 
correlations to a metallic state with singlet pairing correlations. 
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FIGURES 




FIG. 1. Triangular lattice for the dimer model of k-(BEDT-TTF) 2 X. Each lattice site is de- 
noted by a circle and represents the anti-bonding orbital on a dimer of a pair of BEDT-TTF 
molecules. The large arrows show the two possible ladders. 
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FIG. 2. a) The charge gap of the half-filled ladder as a function of U for various values of r. 
b) The spin gap of the half-filled ladder as a function of U for various values of r. 
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FIG. 3. Schematic phase diagram of the "diagonal" ladder (a) and the "horizontal" ladder (b). 




FIG. 4. The rung-rung singlet pairing correlation function D{£) versus I for the half-filled 
ladder for r = 1 and various values of U. 
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